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£ : 

We investigate the causal structure of the Harada-Iguchi-Nakao (HIN)'s exact solution in detail, 
which describes the dynamical formation of naked singularity in the collapse of a regular spherical 
cluster of counterrotating particles. There are three kinds of radial null geodesies in the HIN 
spacetime. One is the regular null geodesies and the other two are the null geodesies which terminate 
at the singularity. The central massless singularity is timelike naked singularity and satisfies the 
strong curvature condition along the null geodesies except for the instant of singularity formation. 
The cluster dynamically asymptotes to the singular static Einstein cluster in which centrifugal force 
, is balanced with gravity. The HIN solution provides an interesting example which demonstrates 

\Q ■ that collisionless particles invoke timelike naked singularity. 

o : 
o 

O ■ I. INTRODUCTION 

"o 

The cosmic censorship hypothesis (CCH) is one of the most important open problems in general relativity since 
$H . it plays an important role in theories of black hole physics The CCH roughly states that singularities forming 

in gravitational collapse must be hidden behind event horizons and hence invisible to outside observers. Many types 
of gravitational collapse have been studied so far in the context of the CCH. Some of them produce globally naked 
singularities, but they cannot immediately be counterexamples to the CCH because the CCH requires suitable matter 
and appropriate initial conditions. 

The well-known example of spacetime which describes the dynamical formation of naked singularity is the Lemaitre- 
Tolman-Bondi (LTB) spacetime ||^] . The LTB spacetimes describe the gravitational collapse of a spherically symmet- 
ric dust cloud and have an explicit form of the entire spacetime metric. It has been proved that the LTB spacetimes 
have ingoing null naked singularity from the generic initial data |p[^10| . However this model is rather simplified because 
pressure is not taken into account. 

Some aspects of the effect of pressure on the naked singularity formation have been investigated. Several important 
results have been obtained. Ori and Piran investigated the spherically symmetric collapse of a perfect fluid numerically 
under the assumption of self-similarity |ll[] , and analytic discussions based on self-similarity followed it IQ . Harada 
also solved numerically the non-self-similar spherically symmetric collapse of a perfect fluid (l3). Their results are 
summarized as that naked singularity can occur in the spherically symmetric collapse of a perfect fluid for a sufficiently 
soft equation of state. In contrast to isotropic pressure, to include only tangential pressure is turned out to be more 
tractable |l4|,[ll|. Among the solutions with vanishing radial pressure, there is a system of a spherical cloud of 
counterrotating particles, in which the physical origin of tangential pressure is clear |lq-|l9[|. The each particle in the 
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cluster has its angular momentum so that the average effect of all particles is a non-vanishing tangential pressure. 
This model is particularly interesting from a physical point of view because it gives insights into rotational effects on 
the gravitational collapse, without raising terrible difficulties. Harada, Iguchi and Nakao (H1N) have recently analyzed 
singularity occurrence in this model |2(] . |2l)| . Using the mass-area coordinates, which were first introduced by Ori |2^| 
and were applied in this system by Magli p3|] , HIN found a new exact solution that describes dynamical formation 
of massless naked singularity. The results show that the counter-rotation undresses the covered singularity, and in 
particular tangential pressure and rotation may induce the formation of naked singularity. Since the HIN solution is 
given in the mass-area coordinates, the motion of each shell and the global properties of the spacetime are not trivial. 
We study in this paper the HIN solution in detail to make them clear. 



The plan of this paper is as follows. In the next section, we review the HIN solution. In Sec. [II, we study null 
geodesies in the HIN spacetime which arc necessary to determine the causal structure of spacetime. It is discussed in 
Sec. IV that the central singularity is timelike. In Sec. |v|, we will see that the collapse asymptotes to the singular 



static Einstein cluster. Section VI is devoted to conclusions. We use units with c = G = 1 and follow the sign 



conventions of the textbook by Misner, Thorne and Wheeler about the metric, Riemann and Einstein tensors [|24| . 

II. THE HIN SPACETIME 
A. The HIN solution 

The HIN solution is the spherical cloud of counterrotating particles which is marginally bound. The specific angular 
momentum L(r) of each particle at comoving radius r equals to 4F(r) pp| , where F(r) is the conserved Misner-Sharp 
mass function pa]. Using comoving coordinates, the line element for this system is reduced to 



ds 2 



, 2 ^dt 2 + R' 2 (t 7 r) (l + ^^l)dr 2 + R 2 (t,r)(de 2 + S in 2 ed^), (1) 
V R 2 (t,r) J 



where v(t,r) and R(t,r) satisfy the following set of coupled partial derivative equations: 

16F 2 

R { R^l^) R '^ < 2 > 

V = e-»R(t,r)=- . (3 ) 
y/R(R 2 + 16F 2 ) 

The prime and overdot denote the partial derivatives with respect to r and t, respectively. The energy density 
e(t, r) = — Tj and the tangential pressure H(t, r) = T® = are given by 

e ^ = lJ^ (4) 

Using the coordinate transformation to — F(r) and R = R(t, r), the HIN solution is given in the mass-area coordinates 
(to, R) by 

ds 2 = -Adm 2 - 2BdRdm - CdR 2 + R 2 {d6 2 + sin 2 6dc/) 2 ), (6) 
where A, B and C are given by 



B \R-4m\\l 2m 1 (8) 
= J_ = R(R 2 + 16to 2 ) 

V 2 2m(R - Am) 2 ' 1 ' 

1\2 



\J H(m, R) = VB 2 - AC 



(F- 1 ), to ((F- 1 ) 2 + 16to 2 ) 
— Am\v / 2mF~ 1 
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R 2 - l&mR + 144m 2 R ^, y/R + 2^E 
=— - — \— +4V21n- 

3V2m(i?; - 4m) V m 



4m) 

- lGmF- 1 + 144m 2 I F^ 1 



3V2m(F- 1 - 4m) 




4^2 In 



(10) 



We denote the inverse function of F(r) as F 1 = F 1 (m) and have set R = F 1 (m) on the initial spacelike 
hypersurface, which corresponds to i?(0, r) — r. The energy density in the mass-area coordinates is given by 



e(R, m) 



VR 2 + 16m 2 



AnR 3 \V\VH 



(11) 



We can assume that the metric functions are C°° class with respect to the local Cartesian coordinates at least in the 
neighborhood of the center r = before encountering a central singularity. By this assumption, the metric variables 
in the comoving coordinates are expanded near the center as 



v(t, r) = I* (t) + v 2 (t)r 2 + v 4 (t)r 4 + 0(r 6 ), 
R{t, r) = R x {t)r + R 3 (t)r 3 + R 5 {t)r 5 + 0(r 7 ). 

Then from Eq. ([|), the arbitrary mass function F(r) is also expanded as 

F(r) = F 3 r 3 + F 5 r 5 + F 7 r 7 + 0(r 9 ), 

and the inverse function _F _1 (m) is approximately given by 



^ _1 (m) = r = 



1/3 



AF 5 2 - 3F 3 F 7 
9F n/ 3 



m 5 / 3 + 0(m J / 3 ), 



(12) 
(13) 



(14) 



(15) 



where we expand F _1 (m) up to the order of m 5//3 because all terms up to this order are needed for later calculations. 
We can set va(t) = by using the rescaling freedom of the time coordinate, and Eq. (||) yields v%(t) — . Thus the 
leading order of v(t, r) is given by 



For R(t,r), we obtain 



where to is given by 



v{t,r) = 



R(t,r) 



4i V 4 , 



Ri\t) 



r 4 + 0(r b ). 



t 



2/3 



r + 0(r 3 ), 



to = TT\ -FT 



1 / 2 
3V F~z 



(16) 



(17) 



(18) 



The break down of these expansions means the appearance of central singularity. Since the break down is given by 
Ri(t) = 0, the singularity formation time is t = to- 

On the initial spacelike hypersurface, the regularity requires R(t — 0,r) = r > AF = 0(r 3 ) in a sufficiently small, 
but finite region around r — 0. From Eq. (|^), R(t,r) is a monotonically decreasing function with respect to t, and 
then each initially collapsing shell approaches R = 4F. Using proper time r(t, r) = J e v dt of each shell, the behavior 
of this approach is 



R — AF cx exp i 



8FJ 



(19) 



This behavior shows that R approaches 4F asymptotically and thus R is always R > 4F if it initially holds. Since 
the trapped region is given by < R < 2F, the region around the center is not trapped eternally. Hence the central 
singularity is globally naked. 
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B. Algebraic root equation method 



The nakedness of the central singularity in the HIN spacetime is also confirmed by examining the algebraic root 
equation, which probes the existence of outgoing radial null geodesies from the center |§|- |l0| , p3|| . To obtain the root 
equation, consider the radial null geodesic equation in mass-area coordinates 



dR 
dm 



c 



J T y/H(m,R), 



where Jzp is given by 



\V\ 



R 



VR 2 + 16F 2 



T\V\ = 



O 2 



l)V2z 



(20) 



(21) 



and we have defined z = R/4m(> 1) for later convenience. The upper and lower signs refer to outgoing and ingoing 
null geodesies in the collapse phase, respectively. Hereafter we use this sign convention. To investigate the behavior 
of null geodesies near the center, we define 



x(m) 



R 



2m al 

where a is a constant. Applying l'Hospital theorem to Eq. (|2^), we obtain the root equation 



(22) 



x 



lim 



+o 2a dm 



lim 



id 



3 ( 1 -°)/ 2 (i -2m 1 - a ) 



^/H(m,2x to") \xl /2 T m {1 ~ a)/2 (x - 2m 1 - 01 ) 



(23) 



where xo = x(0) is introduced. If we find a consistent set of xo and a, we have null geodesies which behave as Eq. 
( f22| ) terminating at the center to = 0. The root equation method picks up only the geodesies behaving as Eq. (g^). 
To find the another possible null geodesies, we must solve the null geodesic equation. 



III. NULL GEODESICS IN THE HIN SPACETIME 

From the detailed analysis of Eq. (p3]) , we can prove that all possible values of a for positive finite xq are 1/3, 7/9 
and 1. In reality, these three values of a correspond to three different types of null geodesies in the HIN spacetime, 
and it will be proved in this section that there are no other null geodesies in the HIN spacetime. We will see below 
the derivations and properties of these three types of null geodesies in detail. 



A. a = 1/3 : regular null geodesies 



The null geodesies with a = 1/3 correspond to regular null geodesies. For this case, the value of xq for a = 1/3 
cannot be determined by the root equation method. We will see that the regular null geodesies actually have a = 1/3. 
From Eqs. (|l|), (|l6| ) and (17), the radial null geodesic equation for lowest order is 



dt ^ 
dr 



1 

to 



2/3 



(24) 



Inserting the solution t = t(r) of Eq. ( p4[ ) to Eq. (|l7j), we obtain the regular null geodesies in the mass-area coordinates 



t(0) 



R{m) M 1 _ 



to 



2/3 



1/3 



(25) 



where t(0) < to is the time when the null geodesies arrive at or emanate from the regular center r = 0. The 
arbitrariness of Xq for a = 1/3 comes from the arbitrary constant t(0) of Eq. (p5h. 
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We can find that the scalar curvature is finite at the center along these null geodesies (see Appendix [A]) . In fact, we 
can prove that all the null geodesies with a = 1/3 terminate at or emanate from the regular center. The integration 
of Eq. (||) by using the proper time r(t, r) is 

r(t, r) = c(m) - 4m f xMfJ^Uds, (26) 

where Ro{m) is the initial area radius i?(0,r) = F~ l (m) and c(m) is an arbitrary function of m. For the regular 
spacetime, v(t,0) — vo(t) = implies that we can set r = t on m = 0. Then the singularity appears onm = when 
the proper time is r = to- The elliptic integral of Eq. ( |26] ) is explicitly integrated by approximating the integrand. 
Ro(m)/4m 1 is satisfied around the regular center m — and if we consider the condition z>l, the integrand of 
Eq. ( p6| ) is approximately y/2s. Then, r for each shell becomes 

t ~™ + *-t(i£>T- (27) 

We consider the null geodesies of Eq. ( g^ ) for a = 1/3. From Eq. (p7|), the proper time along these geodesies is 

T = c(m) + t -^x 3/2 . (28) 

For the regular shell motion of Eq. ( |l7| ) , Eq. ( f27| ) becomes 

T = c(m)+t. (29) 

Taking into account t = < on m = 0, we have c(0) = and thus Eq. ( ^8| ) implies r < to at m = 0. Therefore all the 
null geodesies behaving as R cx m 1 / 3 terminate at m — before the singularity appears. 



B. a — 7/9 : the earliest singular null geodesic 

The sing ular null geodesic with a — 7/9 for both ingoing and outgoing geodesies was found by Harada, Iguchi and 
Nakao EQ] . The behavior of the null geodesic as 

R ps 2x m 7/9 , (30) 
is not regular and thus it terminates at the central singularity. The coefficient Xq is given by 



— l^^)'' (31) 



/ 24F 3 2 _F 5 



for F5 < 24F3 2 . Note that F5 < 24F 3 2 is the same as the requirement of no shell-crossing singularity and this 



condition holds if e(0,r) is a non-increasing function of r. It is found that Eq. (27) reduces to r = to at m — for 
■| < a < 1. Thus all the null geodesies of | < a < 1, particularly a = 7/9, terminate at the singularity. The scalar 
curvature diverges along the null geodesic with a = 7/9 (see Appendix Ah. 

We will see below that there is only one ingoing or outgoing null geodesic with a — 7/9 and that this null geodesic 
is the first one which arrives at or emanates from the singularity at t = to- The null geodesic equation in mass- 
area coordinates has a singular point at (m, R) — (0,0). To make the singular point tractable, we introduce new 
coordinates as 



X = m^ 9 , (32) 
R_* 

,7/9 



R \ 3/2 

7-770 , (33) 



and then Eq. (pQ) becomes 



+ £(tf-A) = A<Mx,0), (34) 
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where is given by 
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(*, = 2 (3X 2 ^ 1/3 Vh(xT)Jt (X. *) " <?) , (36) 

and we have introduced a parameter < A < oo. 

The form of Eq. (|5|) is similar to the form of the null geodesic equation of the LTB spacetime given in |j]|7]] . One 
may now follow to Christodoulou's argument to show the existence and uniqueness of a continuous solution of Eq. 
(p4|) H . It is sufficient to consider the null geodesic equation in the neighborhood of the center \ = since we are 
interested in the radial null geodesies terminating at x = 0. There is no singular point at \ > where $ is strictly 
positive by definition. Problems appear when we consider the center \ = 0. Wc expand Jz^ and vH around x = 
using Eq. (|l|), 



and from this expansion, ip^ is 

MX,#) « + ViWx + V^(^)x 2 + 0(x 3 ), (39) 
where the coefficients of each order are 

Vx) - l2V2x 3 /2 , 



V2 V * 
^(^ = -12^/3 ( 6 + ^ 

If we choose the parameter A to be A = Ao = (2xo) 3 / 2 , ^ T is at least C 1 in x > 0, d > 0. We can apply the contraction 
mapping principle to Eq. (|34|) to find that there exists the solution satisfying $(0) = Ao, and moreover that it is 
the unique solution to Eq. (|34|) which is continuous at x = 0. This solution with $(0) = Ao exactly agrees with the 



geodesies of Eq. (30). The proof is presented in Appendix H. Therefore there is no other solution with < 1?(0) < oo. 



In other words, another possible solution must be $(0) = or oo. 

We consider possible solutions with i?(0) = oo. When i?(0) = oo, Eq. ( |34| ) is approximated around x — by using 
Eqs. ( p7| ) and ( p8[ ) which are valid even in this limit as follows: 

d ± « -2* (41) 

The integration of this equation gives $(x) °c 1/x 6 - This behavior coincides with the regular null geodesies of Eq. 
( p5| ) because all the geodesies behaving as R oc to 1 / 3 terminate at the regular center as we have shown before. 

It is important that R(t, r) is a monotonically decreasing function with respect to t, and thus that "d decreases as 
t increases when x is fixed. Because of this time direction in the mass-area coordinates and the fact that there are 
no other null geodesies with (2cco) 3//2 < i?(0) < oo, the geodesic with $(0) = (2cco) 3 ^ 2 is the first null geodesic which 
arrives at or emanates from the appeared singularity. Hence we conclude that the arrival or emanational time in 
comoving coordinates is the singularity formation time t = to- 

C. a = 1 : later singular null geodesies 

There are also null geodesies with a = 1. x(m) is expressed for a = 1 by a non-analytic function 

x(m) w 2 + 2exp ( — ■= ] , (42) 



G 



where D is the positive constant which parameterizes the null geodesies. We find that the scalar curvature diverges 
at the center along these null geodesies as is seen in Appendix [A]. We will see these null geodesies in detail below. 

In the coordinate these null geodesies are described by solutions with $(0) = if they exist. We first search the 
solution with z(0) = oo and $(0) = 0. Under these conditions, Eq. (|2(]) reduces around to — to 



dm Am 

For z 3 / 2 <C m^ 1 / 3 , we can neglect the second term of Eq. ([y|) and immediately integrate it. However, the solution 
contradicts z 3 / 2 <C m -1 / 3 because it is given by z oc to -2 / 9 . For z 3 / 2 oc to -1 / 3 or z 3 / 2 3> to" 1 / 3 , there are consistent 
solutions. The solution of Eq. (f43| ) for z 3 / 2 > m" 1 / 3 is 

' (44) 



m 2 / 3 ' 

which corresponds regular null geodesies with a = 1/3. When z 3 / 2 oc to -1 / 3 , the integration gives 

z oc to~ 2/9 , (45) 

which corresponds to the unique null geodesic with a = 7/9. These two solutions satisfy the condition z(0) = oo, but 
do not satisfy $(0) = 0. Consequently there is no solution which satisfies these two conditions. 

We consider the geodesies with z(0) < oo and i9(0) = 0. Using 1' Hospital theorem, we obtain a restriction on 

JzfVH as 

dR 

4z(0) = lim — = lim J^Vh < oo. (46) 

m^O dm m->0 T 



Though Jzp is strictly positive at m — as long as z(0) > 1, \/ H diverges to positive infinity at m = as 

VH = 4V2 ( -4 + 7 ^ + 0(m) ) + l^ ^ 2 - 4 - + 9 ) + 4V2 In 4±i, (47) 



where we have defined 



^^4 /2 , (48) 
1440F 3 4 + 48 F 3 2 F 5 - 17 F 5 2 + 12F 3 F 7 



7= — ° 2 - /fi - — . (49) 

216F 3 23/6 



Thus, there exists no solution with the boundary condition 1 < z(0) < oo. However it may be possible for null 
geodesies to satisfy Eq. ( |4^ ) only if z(0) = 1. 

We introduce the new coordinate y = z — 1 to study the solution y{m) with y(0) = 0. By expanding Eq. j20| ) 
around y — 0, we obtain the consistent solution with y(0) = in the lowest order approximation. We give here the 
null geodesic equation which is expanded up to the enough order so that the difference between ingoing and outgoing 
null geodesies appears in the expansion; 

^L^JL(-lny+4-+5TV2 + 1 m^) . (50) 
dm 2m \ to 1 '' 3 / 



It can be integrated to 



y w exp 



^ + ^ + S T V2 + 3 -lm^), (51) 

TO TO ' 5 / 



where D > is an integration constant and S is defined as S = 2 In 2 — | . We will see in Sec. |v| that the constant D is 
related to the time when the null geodesies terminate at the center to = 0. The geodesies of Eq. ( |5l| ) satisfy y(0) = 
and they are classified into the geodesies with a = 1. To obtain the proper time when the null geodesies terminate at 
the singularity, we consider Eq. (|26| ) again. By the approximation of the integrand, Eq. (^6|) reduces around to = 
to 

t ^ c{m) - 4m(^2\ny{m) - . (52) 

Then the proper time when the null geodesies terminate at to = is r = t and it means that the null geodesies 
terminate at the central singularity. As a result, it is concluded that all the solutions with z?(0) = are the solutions 
which terminate at the singularity with y(0) = 0. 
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IV. CAUSAL STRUCTURE OF THE HIN SPACETIME 



A. Timelike singularity 

It is important that Eq. (|5l| ) includes all the null geodesies emanating from the center after the singularity appears. 
We construct double null coordinates by using Eq. (^Tj) and study the causal structure of the spacetime which is 
covered by these null geodesies. From Eq. (|5l|), we introduce double null coordinates (u,v) which satisfy around 
m = 

2V2u^ ^i(-\ny+^- + S-V2 + ^m 1 A 7 (53) 
\ m i / J 5 / 

2V2«« ^T l (-\ny+^- + S+V2 + ^m 1 A . (54) 



Then, m and y are expressed in the null coordinates by 

\fm w v — u, (55) 

y « exp ^ '- + + S + -± (v - ufl* . (56) 

I v — u [v — u) z ' a 5 ; 



We are considering a sufficiently small, but finite region around m = 0. Then Eq. (55) restricts u and v to u rs v. 
From Eqs. (pi) and (pq), dm and dR are given by 



dm ss 2(w — u)(dv — du), (57) 
dR-Adm^ -4V2y (^v + u - s/2(3(v - u) 1/12 + (y/2 + l)(v - u) + ^-"f (v - u) 5/12 ^j du 

+AV2y (^v + u- V2f3(v - u) 1/12 + (y/2 - l){v - u) + ^7 (v - u) 5/12 ^ dv. (58) 

Inserting Eqs. (p7| ) and ( p8[ ) into Eq. (||) and expanding the metric functions A, B and C around (m,y) — (0,0), we 
obtain the line element in the double null coordinates as 

ds 2 w -512(v - ufdvdu + R 2 {d8 2 + sin 2 8d(j) 2 ). (59) 

The details of the calculation are given in Appendix^. In the null coordinates, the central singularity is represented 
at u — v which corresponds to the center m — because of Eq. (j55|). It is found that the world line of u = v is 
timelike. Hence the central singularity in the spacetime which is covered by the null geodesies given by Eq. ( pl| ) is 
timelike. 



B. Penrose diagram 



Here we summarize the obtained results. The whole of the spacetime is covered by three types of null geodesies 
which can be classified by the value of a into a = 1/3, 7/9 and 1. The null geodesies with a = 1/3 emanate from 
or terminate at the regular center and they are parametrized by one parameter. There is only one null geodesic with 
a = 7/9 and it is the earliest one which emanates from or terminates at the naked singularity. The null geodesies 
with a = 1 emanate from or terminate at the timelike singularity and they are parametrized by one parameter. From 
these results, it is now possible to draw the conformal diagram of the HIN spacetime (see Fig. |l|). For comparison, 
we also present the conformal diagram of the LTB spacetime with naked singularity (see Fig. ||). It is found that the 
effect of counterrotation makes the singularity timelike in this model. 



C. Curvature strength of the singularity 

We should note the curvature strength of the singularity in the HIN spacetime. According to Tipler and Krolak, we 
classify the curvature strength whether the singularity satisfies the strong curvature condition (SCC) or the limiting 
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focusing condition (LFC) 26 2^]. Harada, Nakao and Iguchi studied the SCC and the LFC for spherically symmetric 
spacetimes with vanishing radial pressure p8| . Applying theorems 1, 2 and 3 of their paper to the HIN spacetime, we 
can know the curvature strength along each null geodesic. For a — 7/9, not the SCC but only the LFC is satisfied 
along the null geodesies. For a = 1, the SCC is satisfied along the null geodesies. In conclusion, the naked singularity 
is relatively weak at the formation but becomes strong after that. 



V. ASYMPTOTIC BEHAVIOR 



It is worth while examining the asymptotic behavior of the spacetime. We concern with the null geodesic behavior 
in the comoving coordinates corresponding to Eq. (pl|). To obtain an insight, we consider the asymptotic behavior 
of the HIN solution in comoving coordinates. Since the time coordinate t at the center does not proceed after the 
singularity formation, we must introduce new time coordinate which no longer agrees with the proper time at r = 0. 
We denote the new time coordinate as T and consider that the time t in Eqs. (|j), (||) and (|J) is replaced by T. 

We take the limit R = AF because R{T, r) approaches AF(r) asymptotically. From Eq. the metric in this limit 
is given by 

ds 2 = if dT 2 + 2{AF\r)) 2 dr 2 + (AF(r)) 2 (d0 2 + sin 2 Odcf 2 ). (60) 

K(T) 

The function K(T) is an arbitrary function of T which comes from the integration of Eq. (Q). For simplicity we set 
K(T) = K = constant by rescaling the time coordinate. The energy density in this limit is 



1 



256nF 2 ' 



(61) 



This solution is a static system of counterrotating particles, which is called the Einstein cluster. In particular, the 
solution has timelike naked singularity at the center. 

To study the asymptotic behavior of the shell motion for fixed r, we set perturbed quantities e(T, r) and d(T, r) for 
the coordinates (T, r) around R = AF as, 

R(T,r)=AF(r)[l + e(T,r)}, (62) 

e 2,(T,,) = ed (T,r) BIA « [1 + e(Tj r) + d(T , r)] . (63 ) 

Inserting these quantities into Eqs. (||) and (|J), we obtain following perturbed equations in linear order, 

e(T,r)=--4M=, (64) 



A^KF^r)' 

F'(r) 
F(r) 



F'(r) 

d'(T,r) = -— ff e(T,r). (65) 



e(T, r) = E(r) exp , (66) 




The solutions of these equations are 



d(T,r) = - J e (T,r)^dr + G(T), (67) 

where E(r) and G(T) are arbitrary functions of the comoving radius r and time T respectively. We have implicitly 
assumed that the time T is very large compared to the singularity formation time. However, Eqs. ( |66| ) and ( |67| ) imply 
that this perturbation scheme is valid soon after the singularity formation time as long as we consider the region in 
which the radius r is very small. 

In the asymptotic region, the proper time becomes 

T (T,r)= 2\ — [T + const. + O (e(T, r)) ] . (68) 
V K 
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It is found that the asymptotic behavior of Eqs. (|66|) and (|68|) completely coincides with the already obtained behavior 
of Eq. (19). Using these results of perturbation, the null geodesic equation in the comoving coordinates is 



— = ±\—(2F' + 2Fe' - F'd 
dr V r \ 



I2K 



and it is immediately integrated to 



T(r) = T(0) ± V2K (iy/W) + J ^= [2*V ~ F'd] 



(69) 



(70) 



Integration constant T(0) is the time when these null geodesies terminate at or emanate from the singularity, and 
this parameterization of the null geodesic family corresponds to the fact that the singularity is timclike. From Eqs. 
(p2fi, (pq) and (fTO), we obtain the null geodesies in the mass-area coordinates corresponding to Eq. (179), 



y= — -l = E(r)exp( ^£L 



1 



dr 



= [2Fe'-F'd] 



2^2F{r)J y/FfSj' 

Comparing this result to Eq. (El]), the parameter D relates to the arrival time T(0) as T(0) = 4Dy/K. 



(71) 



VI. CONCLUSIONS 



We have studied the causal structure of the HIN spacetime and it was shown that the central massless singularity 
of this spacetime is timelike. To show this fact, we have investigated the radial null geodesies in detail. The null 
geodesies are classified to three types, a = 1/3, 7/9 and 1, by their dependence of R on m near the center. One 
is regular and the other two are singular. The classification of the singular geodesies corresponds to their arrival or 
emanational time at the central singularity. The a — 7/9 type is the earliest singular null geodesic which arrives at 
or emanates from the singularity at its formation time to, while the geodesies with a — 1 arrive at or emanate from 
the singularity after its appearance. 

We have shown that singular null geodesies with a = 1 exactly parametrized by their arrival or emanational time 
and that there is only one set of ingoing and outgoing geodesies for each parameter. This fact shows that the central 
singularity has timelike property. We have also constructed double null coordinates around the central singularity 
from the null geodesies with a = 1. The line element in this double null coordinates shows that there is timelike 
singularity in this spacetime. We have considered the asymptotic behavior of this spacetime after the singularity 
appeared in comoving coordinates. It gives us understanding of the null geodesic behavior in mass-area coordinates, 
and of the parametrization of the null geodesic family. The curvature strength of this singularity was also investigated. 
The LFC is satisfied along the null geodesies with a — 7/9 and the SCC is satisfied for a = 1. The curvature strength 
of the naked singularity is relatively weak at the formation and it becomes strong after that. 

In summary, the HIN solution describes a dynamical formation of timelike naked singularity, that is, the birth of 
timelike singularity. The solution dynamically tends to the static singular Einstein cluster as time proceeds. Though 
the HIN system is simply composed of collisionless particles, the collapse leads to the nontrivial causal structure. 
It implies that the effects of rotation and tangential pressure play important roles in the final stage of collapse, 
particularly the singularity formation. 
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APPENDIX A: SCALAR CURVATURE ALONG THE NULL GEODESICS 



Singularities are boundary points of spacetime where the normal differentiability breaks down. If the energy density 
or the curvature invariant diverges at boundary points, the points are singularities. In this appendix, we give the 
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scalar curvature in the HIN spacetime along the radial null geodesies terminating at the center r = 0, and show 
that the center is singular along the geodesies with a — 7/9 and 1. 
The scalar curvature i?^ of the HIN spacetime is given by 



Rf, 



8-kR 2 



V2 



" i? 2 + 16F 2 (R-4m)VmRVH' 
Along the geodesies R — 2xm 1 / 3 , which are shown to be regular null geodesies, R^ is given by 

3 



(Al) 



(A2) 



and it is finite at m = 0. On the other hand, when we consider R£ along the geodesies, R = 2xom 7 / 9 and R — 4m(l+y) 
of Eq. (51), the scalar curvatures are given by 



Rf, 



28x 3 m 4 / 3 



0(m 



-10/SK 



(A3) 



and 



Rf, 



1 



+ 0(m- 5/2 y) 



i* 64m 2 

respectively. The scalar curvature diverges as m — ► in these cases, and therefore the center is definitely singular 



(A4) 



APPENDIX B: CONTRACTION MAPPING PRINCIPLE 

We prove the existence of the unique continuous solution of Eq. (|34|). It is based on the discussion given by 
Christodoulou for the LTB model M. Consider the differential equation obtained from Eq. ( [m]) by replacing -d in 
^t(Xj^) by a given continuous function -d > 0: 

£E: + £(0-Ao) = AoMx,0(x))- ( B1 ) 

dx X 

Since we have chosen A = Ao, ^^(x> v) is at least C 1 in the strip x G [0, Xi] an d $ G (0, Mi], where we assume xi 



is sufficiently small. The continuous solution of Eq. (Bl) is only the solution with i?(0) = Ao(> 0). This solution is 
given by 

0(X) = X n( 1 + *l S ^T ( s *> ^( S X)) rfs) • (B2) 

We consider the nonlinear map T\ defined as i? = T\ (d). Since we are considering the function $ with $(0) > 0, 
we obtain l\ o (#)(0) = A from expansions of Eqs. ([57) and (pSj). Thus the possible fixed point of T Ao is only one 
which satisfies i?fp(0) = A . To prove the existence of the fixed point, we consider the set consisting of all $ such 
that for x G [0,X2], 

ft- < Hx) < ft+, (B3) 

where the lower and the upper bound are defined by /i_ = Ao — fi2 and /i+ = Ao + /i2 for a sufficiently small /12 . From 
these bounds, the nonlinear map is also restricted in 

ft- <TxM(x) <M+, (B4) 

if we choose a x sucn that 

X < X3 (B5) 

where 
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Ai = sup sup 

o<x<X2 n_<a<fj, + 



(B6) 



The map T\ sends into itself for all \a < m i n {X2,X3j-- Let 



A2 = sup sup 

0<X<X2 ^i_<i3</j 4 

Then we obtain from Eq. (fB2]) for $1, $2 € V^, 



|r Ao (^)-r Ao (^)||<^^||^_^|| 



(B7) 



(B8) 



where || || denotes the supremum norm. If we choose \i < 7/(AoA2), the map T\ g is contractive in V^. Hence T\ 
has a unique fixed point $fp € which is given by 



(B9) 



Therefore we have the unique continuous solution with i?(0) = Xq. 



APPENDIX C: DERIVATION OF EQ. (59) 



In this appendix, we construct double null coordinates by using Eq. (^lj). We define double null coordinates (u, v) 
as Eqs. (p3) and (p#. From this definition, dm and dR are given by 



dm = p u du + p v dv, 

dR = Adm + (q u du + q v dv), 



where 



Pu 
(Ju 



-p v w 2(u - v), 



-4V2y L) + u- V2/3(v - u) 1/12 + (V2 + l)(v - u) + ^7 (v - u) 5/12 J du, 

q v m 4\/2y + it - V2(3(v - w) 1/12 + (\/2 - - u) + ^7 (w - u) 5 / 12 ^J dv. 

The line element of Eq. ([!]) is rewritten by the coordinate transformation from (m, i?) to (it, v) 

ds 2 = Udu 2 + Vdv 2 + Wdudv + R 2 dfl 2 , 

where the metric functions U, V and W are given by 

U = -(A + 8B + 16C)p 2 u - Cq 2 u - 2(B + AC) Pu q u , 
V = -(A + 8B + 16C)p 2 v - Cq 2 v - 2{B + AC) Pv q v , 
W = -2(A + 8B + l6C) PuPv - 2Cq u q v - 2{B + 4C)( Pu q v + p v q u ). 



From Eq. (C3), we obtain 



and then Eq. (|Cq ) becomes 
ds 2 w - 

If the condition 



-U-V-C(q u +q v ) 2 , 
-C(q u + q v ) 2 dudv + U(du 2 - dudv) + V(dv 2 - dudv) + R 2 dVt 2 . 



(CI) 
(C2) 



(C3) 
(C4) 

(C5) 
(C6) 

(C7) 

(C8) 
(C9) 
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C(q u + q v ) 2 ^U,V 



(CIO) 



is satisfied, the double null coordinates would be constructed approximately. 
We expand \f~H of Eq. (Et]) around y = 0, 



where h(m) is defined as 



// = h(m) + V2 f ~ - 41ny + 4 (3 + S) + y + ^y 2 + 0(y 3 ) ) . 



h(m) = 4V2 f-JL + 7m 1 / 3 + 0{m)\ 



Inserting Eq. (Cll) into A, B and C, metric functions U, V and W are calculated. As a result we obtain 



(Cll) 



(C12) 



lim C(q u + q v ) 2 = -512, 

m— >0 m 

u 

lim — = 0. 

m^O m 

V 

lim — = 0, 

m— >o m 



and it satisfies the condition of Eq. (CIO). Hence the line clement around the center is given by 

ds 2 w -512(« - u) 2 dudv + R 2 dVL 2 . 



(C13) 



(C14) 
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FIG. 1. The conformal diagram of the HIN spacetime. A timelike naked singularity emerges at the center at t = to- 
There are one parameter family of both ingoing and outgoing null geodesies which terminate at and emanate from the naked 
singularity. 




FIG. 2. The conformal diagram of the LTB spacetime with naked singularity. There are one parameter family of outgoing 
null geodesies which emanate from the naked singularity but only one ingoing null geodesic which terminates at the naked 
singularity. 
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